AFFINE SYNTHESIS AND COEFFICIENT NORMS FOR LEBESGUE, 
HARDY AND SOBOLEV SPACES 



H.-Q. BUI AND R. S. LAUGESEN 

ABSTRACT. The affine synthesis operator Sc — X^>o X^fegz^ c i,fcV'j,fc ls shown 
to map the mixed-norm sequence space t L (t p ) surjectively onto L p (R d ), 1 < 
p < oo, under mild conditions on the synthesizer ip £ I/ p (R d ) (say, having a 
radially decreasing L 1 majorant near infinity) and assuming J Rd ip dx = 1. Here 
4>j,k(x) = | det aj\ x ^ p il}(ajX — k), for some dilation matrices a,j that expand. 

Hence the standard norm on / 6 L p (R d ) is equivalent to the minimal coeffi- 
cient norm of realizations of / in terms of the affine system: 

ii/ii, « m J \^n 1/p = / = ££ c ^> k \ ■ 

We further show the synthesis operator maps a discrete Hardy space onto 
1 (R d ), which yields a norn 
Riesz kernel sequence {ze}: 



H 1 (M d ), which yields a norm equivalence involving convolution with a discrete 



ll/llffi ~ mf < E E + I E z Z c i*-t\) ■ f = E E c i< k ^ k \ ■ 

(j>o ket d t£i, d j>o kez d ) 

Coefficient norm equivalences are established also for the Sobolev spaces 
W m ' p (R d ), by applying difference operators to the coefficient sequence Cj t k- 



1. Introduction 

This paper studies mapping properties of the affine synthesis operator 

c = {cj, fc } >-» E c i^h k = ^ 
i>o k& d 

where 

tpj,k(x) = | detajl^^ajX - x G M d , 

assuming J Rd ipdx = 1. Affine synthesis arises naturally in harmonic analysis and 
approximation theory, as a discretization of convolution. 
We first explain our notation, and then our goals. 

• The dimension d G N is fixed throughout the paper, as is the exponent p. 
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• The dilation matrices aj are invertible d X d real matrices that are expand- 
ing: 

|| — > as j — > oo. 

(Here || • || denotes the norm of a matrix as an operator from the column 
vector space M. d to itself.) For example, one could take aj = 2 J I. 

• The translation matrix b is an invertible d x d real matrix, for example the 
identity matrix. 

Our goal is to synthesize surjectively onto the classic function spaces of analysis, 
while assuming as little as possible about the synthesizer tp. This will demonstrate 
that the ability to decompose arbitrary functions into linear combinations of the 
translates and dilates ipjk does not require any special properties of ip, even though 
the efficiency of decomposition naturally does depend on such properties. 

Lebesgue space. In |5] we showed 

S : {finite sequences} — * LP has dense range 

for 1 < p < oo, assuming only that ip has periodization locally in LP (mean- 
ing Y^k& d \i>{x — bk)\ G Lf oc ) and has nonzero integral (f Rd ipdx / 0). This 
density of the range of S means that the small-scale affine system {tpj k : j > 
0, k G Z d } spans LP . Notice Strang-Fix conditions are not imposed: the translates 
{ip(- — bk)} k&Z d need not form a partition of unity. And recall from Q that the 
periodization assumption on ip is easily satisfied, holding for example if tp G L p 
has a radially decreasing L 1 majorant near infinity. 

Two natural questions arise from this L p density result: what is the right do- 
main for 5? and does S map this domain onto L p ? We answer these questions in 
Section [3l where Theorems Q] and [2] show that 

S : l 1 ^) — > L p is linear, bounded and onto. 

The domain £ l (£ p ) is the mixed-norm space of coefficients satisfying 

E(Em p ) 1/p <°°- 

j>0 k&Z d 

From the surjectivity of S onto LP one deduces a coefficient norm equivalence 
of the form 

II/IIp ~ inf{||c|| £ i (£P) : / = Sc}. 
(Bruna H Theorem 4] earlier proved the case p = 1.) The constants for this 
norm equivalence are evaluated in Corollary [3l and Corollary 0] shows that in fact 
equality holds under suitable conditions. 

Corollary [5]proves a localized norm equivalence, on a domain Q C M. d . 

Hardy space. Section 0] considers the Hardy space H 1 = H 1 (R d ) = {/ G L 1 : 
f * (x/\x\ d+l ) G L 1 }. Again we ask what the domain of the synthesis operator 
should be, and whether S maps this domain onto H . Theorems [7] and [8] answer 
the question by showing 

S : ^(h 1 ) — > H 1 is linear, bounded and onto, 
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provided ip € L has nonzero integral and is somewhat "nice" (being for example 
an LP function with compact support for some p > 1, or a Schwartz function). The 
discrete Hardy space h 1 here is defined in Section 0] by convolution of sequences 
against a "singular" kernel in the ^-variable (which is the discrete analogue of 
convolution against x/\x\ d+1 in the definition of the continuous Hardy space H l ). 

Notice the synthesizer ip cannot belong to H 1 , because it has nonzero integral, 
and so to ensure the linear combination Sc belongs to H 1 we must invoke can- 
cellation properties of the coefficient space h 1 . This approach to synthesis in H 1 
seems natural to us because the analysis operator takes H 1 to h 1 (see below), and 
one would like to reconstruct H 1 using only the data obtained from analysis. 

Corollary [9] deduces a coefficient norm equivalence of the form 

\\f\\ H i PS inf{||c||^ (/l i) : / = Sc}. 

Sobolev space. For synthesis into the Sobolev space W m,p (M. d ), we establish 
boundedness and surjectivity in Theorems ITTI and [121 of Section [5] The sequence 
space from which we synthesize involves difference operators with respect to the 
A>index, which act as discrete analogues of differentiation. 

Discussion. Our results on surjectivity of the synthesis operator seem to be 
qualitatively new — they rely on a method of "scale-averaged convergence" that 
we developed only recently for LP in Q. Boundedness of the synthesis operator 
has of course been studied before. 

We also consider boundedness of the affine analysis operator at each scale j, in 
other words, boundedness of / i— > {(/, (pj,k)}kez d from LP to £ p , from H 1 to h 1 
and from W m,p to w m,p . See Propositions [T6l I2T1 and Irrespectively. We show 
the full analysis operator (over all scales) maps isomorphically onto its range in 
Corollaries [6l [10] and [T33 thus giving coefficient norms in terms of the analysis 
operator. 

If one works at a fixed scale j, rather than considering all scales j > as we do 
in this paper, then synthesis yields a shift invariant subspace of L p . Aldroubi, Sun 
and Tang's p-frame work for such shift invariant spaces [1] is described at the end 
of Section [3] 

Topics we do not pursue in this paper include Gabor systems (modulations and 
translations) and wavepacket decompositions (modulations, translations and dila- 
tions). For some recent work in these areas one can consult lfl2l[T7l[T9ll27l[3Tl . 

The paper is structured as follows. Section|2]establishes notation and definitions. 
Sections [3]-f5]present our synthesis results on Lebesgue, Hardy and Sobolev spaces. 
The proofs are in Sections l8l-IT0"l Appendices treat discrete Hardy spaces, and 
Banach frames. 

Parallel results on Triebel-Lizorkin spaces are discussed in Section [6l Open 
problems when ip has zero integral ( f Rd ip dx = 0) are treated in Section |7J in- 
cluding Meyer's Mexican hat spanning problem for LP and its counterpart for the 
Hardy space. We hope this paper helps contribute towards an eventual resolution 
of these fascinating open problems. 
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2. Further definitions and notation 

1. We use doubly -indexed sequences c = {cj t k}j>o kez. d °^ com pl ex numbers, 
with the norm 

Mew : = \Y1 I c j>i 

when 1 < p < oo. Whenp = oo, define ||c||^i^oo) := ^>o su Pfcez d l c i,fel- Then 

^(F) := {c : ||c||^i^p) < oo} 

is a Banach space. 

2. Write LP = L^IR ') for the class of complex valued functions with finite 
L p -norm. Given ip G LP and € L 9 , with 

i + i = i 

P 9 

by notational convention, we define rescalings 

^j,k{x) = | det aj\ 1 ^ p tfj(ajX — bk), (j)j^{x) = \ det aj\ 1 ^ q (j)(ajX — bk). 

These rescalings preserve the L p -norm HV'j.fcllp = 1 1 ^ lip an d the L 9 -norm H^^fcllg = 
\\(j>\\ q , respectively. Alert: the notation ipj^ conceals its dependence on p. 

3. The synthesis operator is 

Sc = Sj, tb c = c i- fc ^'. fc - 
i>o kei d 

Our theorems will specify acceptable domains for this operator, and will explain 
the sense in which the sums over j and k converge. Occasionally we will synthesize 
at a fixed scale j by writing 



SjS = ^2 s k ip jjk , 



for sequences s = {sk}ke'L d - 

4. The analysis operator at scale j is 

T 3 f = T h<p f = {\detb\{fA h k)} k& d - 

That is, Tj maps a function / to its sequence of sampled (^-averages at scale j. The 
full analysis operator simply combines these sequences as 

r/ = r / = {|det6|(/,^ )fc }} J . >O)fe6Zd . 

Our analysis and synthesis operators depend implicitly on the exponent p, through 
the definitions of cpj^ and tpj^. 

5. The periodization of a function / is 

Pf(x) = | det b\ ^2 f( x ~ hk ) for x G Rd - 
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6. Write C = [0, l) d for the unit cube in R d , and C = (-1/2, l/2) d for the 
centered open unit cube. We regard C as consisting of column vectors and Co as 
consisting of row vectors, as the context will always make clear. 

3. L p results 

First we obtain boundedness of the synthesis operator, when the periodization 
of the synthesizer belongs locally to LP . This was already observed by Aldroubi, 
Sun and Tang [1 , formula (2.3)], but we give a proof in Section [8Tl anyway, to keep 
the paper self-contained. 

Theorem 1 (Synthesis into L p ). Assume 1 < p < oo and tp G L p with P\tp\ G 
Then S : £ l (£ p ) — * LP is bounded. More precisely, if c G ^(iP) then the series 

Sc = Yl Yl c i,ki>j,k (i) 

j>o fcezd 

converges in LP in the sense that 

the sum over k in (Q]) converges pointwise absolutely a.e. to a function in 
LP, and the sum over j converges absolutely in LP, 

and furthermore 

\\Sc\\ p < IdetftriP^HlLP^ollcll^^). (3) 

After proving the theorem in Section [8TTI we also give examples to show i l (i p ) 
is the "correct" domain for S, when the synthesizer tp has nonzero integral. When 
p = 2 and tp has zero integral, we point out that synthesis can be bounded on the 
larger domain I 2 (I 2 ), for wavelet and affine frame generators. 

Remember that Sc in (GQ) depends implicitly on the value of p, through the renor- 
malization factor | det aj\V p in the definition of ifijk- This dependence would be 
problematic if we synthesized into more than one L p -space at a time, but we will 
not. 

The hypothesis that the periodization of | ip\ be locally in LP is rather weak, and 
is easily verified in many cases. For example when p = 1 it holds for all tp G L 1 . 
And for p > 1 it holds when tp G LP has compact support or when tp has a bounded, 
radially decreasing L 1 -majorant, or when tp equals a sum of such functions. On 
the other hand, P\tp\ G L p oc can hold even when tp does not decay at infinity. See 
(5] §3.1] for all these observations. Thus Theorem [TJ improves somewhat on earlier 
boundedness results (which go back as far as 11281 l36l 1371 ) because tp need have 
neither compact support nor decay at infinity. 

In the case p = 2, the periodization hypothesis on tp can be weakened to just tp G 
L 2 with P(\tp\ 2 ) G L°°. Specifically, one has bounded synthesis with \\Sc\\2 < 
||-f(|' i / ? | 2 )||oo 2 || c ll^ 1 (^ 2 ) where this last periodization is taken with respect to the 
lattice T, d b~ l . This estimate is proved in (9l Theorem 7.2.3], by showing that 
P{\tp\ 2 ) G L°° is exactly the condition for the integer translates of tp to form a 
Bessel sequence (that is, to satisfy an upper frame bound). 
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For all p one might ask whether the periodization assumption on \tp\ in Theo- 
remCQcan be weakened to just tp G L 1 n L p . We do not know. 

Next we show the synthesis operator is surjective, when 1 < p < oo. In other 
words, we show every / G L p can be expressed by a series of the form (Q}. 

Theorem 2 (Synthesis onto L p ). Assume 1 < p < oo amc? tp £ LP with P\tp\ £ 

r/je« 5* : — > L p /s o/?ew, ant/ surjective. Indeed, if f £ LP and e > ?/ze« 
a sequence c G £ (£ p ) emfs swc/j f/ja? Sc= f with convergence as in a«<i swc/j 

llcll^^ldetbl^H/llp + e. 

Section [8~4l has the proof. The integral of tp is well defined, in the statement of 
Theorem |21 because the assumption P\ip\ G L p oc implies P\ip\ G Lj oc and hence 
ip G L 1 . 

We are not aware of any prior general work on surjectivity of the synthesis op- 
erator, when tp has nonzero integral. The closest work seems to be Filippov and 
Oswald's construction in EUl Theorems 1 and 3], ED . of "representation systems" 
by which every / G LP can be written as a convergent series Sc = f. This looks 
like surjectivity, but the drawback is that their result yields no control over the size 
of coefficients in the sequence c, and thus it is unclear what the domain of S really 
is. Further, Filippov and Oswald work only with isotropic dilation matrices. 

We also mention the density results coming from Strang-Fix theory (discussed 
in §3], although note Theorem |2] holds without needing the Strang-Fix hy- 
potheses). When tp has zero integral, the wavelet theory (9l [13] |29l [32l and related 
phi-transform theory (described in Section © provide sufficient conditions for ob- 
taining frames, orthonormal bases, and unconditional bases, provided the large 
scales j < are included in the synthesis. These conditions all ensure surjectivity 
of S on suitable domains. But the zero-integral case also raises intriguing open 
problems, discussed in Section|7] 

For the special case p = 2, we remarked above that the condition P(\ip\ 2 ) G L°° 
implies bounded synthesis. We suspect it also implies surjectivity onto L 2 , pro- 
vided 1*01 is continuous near the origin and tp(0) = 1. These conditions certainly 
guarantee the ipj^ span L 2 , by a result of Daubechies [ 13 , Proposition 5.3.2], and 
thus S has dense range. One would like to improve this to full range. 

Remark on non-injectivity. The synthesis operator is not injective, and indeed has 
a very large kernel. For example, we could discard the dilation a\ (in other words, 
discard all terms with j = 1 in the sum defining Sc) and still show S maps onto 
L p , by applying Theorem [2] with the remaining dilations {02,03, • • •}■ 

Equivalence of the L p and norms follows immediately from Theorems Q] 

and ED 

Corollary 3 (Synthesis norm for LP). Assume 1 < p < 00 and ip £ L p with 
P\ip\ G L p oc and J Rd tp dx = 1. Then 

II/IIp ~ inf (l|c||£i(£P) ■ f = Sc as in ©} 
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for all f G L p . Explicitly, 

|det6|||P|^|||^ (6C) ||/|| p < inf {||c||^ P) : / = Sc as in ©} 

< | det^ll/ll,. 

For p = 1, the corollary was proved by Bruna [4, Theorem 4]. His duality 
methods apply without needing our assumption that the translations be restricted 
to a lattice. Recall when p = 1 that the periodization condition P\tp\ G LP, is 
superfluous, holding automatically for all tp G L 1 . 

Remark on norm equivalence. As soon as S maps orcfo L p , we know the map 
S : I 1 (£ p ) I ker S —> LP is a bounded linear bijection. (Here we use the canonical 
norm on the quotient space: inf c / e k cr 5 ll c + c ' lU 1 ^)-) Then the inverse map S^ 1 
is bounded by the closed graph theorem, giving equivalence of the £ l (£ p ) and L p 
norms like in Corollary [3] But Corollary [3] goes further, for it provides an explicit 
upper bound for the norm equivalence, based on explicitly estimating the norm of 
in Theorem 12 

Corollary 0] goes further still, giving actual equality of the norms when tp is 
nonnegative. 

Corollary 4 (Synthesis norm equality). Assume 1 < p < oo and that tp G LP is 
nonnegative. When p = 1 assume J Kd tpdx = 1, and when 1 < p < oo assume 
Pip = 1 ( which implies f Rd tp dx = 1). 
Then for all f € LP, 

II/Hp = |det6r 1/9 inf {l|c||^(#>) ■ f = Sc as in®} . 

The constant periodization condition Pip = 1 in this corollary says that the 
collection {| det b\ip(x — bk) : k G Z> d } of translates of ip is a partition of unity. 
Examples of such tp (when b = L) include the indicator function 1c and convolu- 
tions of this indicator function with any nonnegative function having integral 1. 

We next "localize" Corollary [3] to an open set Q C M. d . To that end, we say a 
sequence c = {cj,k}j>o,k£Z d i s adapted to and tp if spt(tpj t k) C f2 whenever 
Cj fe ^ 0, or in other words if Cj^ = whenever spt(V>j,fc) fl J7 C 7^ 0. The point of 
this definition is to ensure Sc = on the complement of $7. 

Corollary 5 (Synthesis norm for L P (J7)). Assume $7 C M d is open and nonempty, 
take 1 < p < 00, and suppose tp G L p (M. d ) is compactly supported with f Rd tp dx = 
1. Then for all f £ L p (n), 

II/IIlp(c) ~ m f { IMI^ 1 ^) '■ f = Sc as in (f2]), and c is adapted to O and tp} . 

The constants in this norm equivalence are the same as in Corollary [3j and so 
they depend on tp and b but are independent of Q. The corollary is proved in 
Section [83] 

We turn now to the analysis operator, which also yields a coefficient norm. 

Corollary 6 (Analysis norm for LP). Assume 1 < p < 00, and take an analyzer 
(p G L q with P\(p\ G L°° and J Rd 4>dx = 1. 
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Then for all f € LP, 



) 



l/p 



11/11 



P 



\\Tf\ 



£<»(#>) = SUP V \(f,4>j,k)\ P 
j>0 VifeeZ* 



det 6|. 



In other words, the analysis operator is linear, bounded and injective from IP 
onto its range in the mixed norm sequence space £°°(£ p ). To explain the appearance 
of £°°(£ p ) in the corollary, note the analysis and synthesis operators are adjoints, 
with Tj, : L p -> l°°(P>) being the adjoint of S4 : £\£ q ) -> L 9 , at least when 
1 < p < oo. Thus the injectivity of analysis in Corollary [6] is equivalent to the 
surjectivity of synthesis in Theorem [2] But we prove the corollary in Section [831 
anyway, for the sake of concreteness and to handle p = 1. 

We close the section by describing Aldroubi, Sun and Tang's work on p-frames 
HI, which is close in subject matter to this paper but has little direct overlap. They 
study the shift invariant range space Vj = Sj(£ p ) of the synthesis operator at a sin- 
gle scale j (whereas our results combine all scales j > 0). Roughly, they showed 
that the reconstruction formula SjTj = identity on Vj (which says functions in Vj 
can be synthesized from their sampled average values at scale j) holds if and only 
if Vj is closed in L p , if and only if Tj is injective on Vj when cj) = tp, if and only if 
V> satisfies a certain "bracket product" condition. This is all carried out in the mul- 
tiply generated case, with synthesizers ipi,...,ip r . They observe that the range Vj 
need not be closed, for example when ip = l[o,i) — l[l,2) i n one dimension with 
1 < p < oo H] page 7]. 



Our Hardy space results assume the dilation matrices are isotropic and expand- 
ing. To be precise, all the results in this section assume that 



for some "dilation" sequence a = {aj}j>o of nonzero real numbers with \ctj\ —>■ 
oo as j — > oo. 

We will recall the Hardy space H 1 , and then construct a discrete Hardy space 
on which the synthesis operator can act. Then we state our synthesis results. 

Hardy space H . Define the Fourier transform with 2ir in the exponent: 



4. Hardy space results 



cij = aJ 




for row vectors £ G M d . Write C d = T{{d + l)/2)vr-( d+1 )/ 2 and 




for the Riesz kernel, so that the Riesz transform of / G L is 



Rf(x) = (/ * Z)(x) = p.v. f f{x - y)Z(y) dy. 
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Then Rf is finite a.e., and is a measurable vector-valued function of x G M. d . 
Notice 

Rf(0 = 

Recall the Hardy space is 

H 1 = # 1 (R a! ) = {/ G L 1 : Rf G L 1 } 

Functions in the Hardy space have vanishing integral: if / € H 1 then i?/ G L 1 
and so Rf is continuous, which implies 

f(0)= f f(x)dx = and Rf(0) = f Rf(x)dx = 0. (4) 

The Riesz transform commutes with dilations and translations, meaning: R(f(ax— 
xo)) = sign(a)(Rf)(ax— xq) when a G M\{0},xo G M. d . But dilation invariance 
fails when a is an arbitrary matrix, which is why we restrict to isotropic dilations 
in this section. 

All these facts about Riesz transforms and H 1 can be found in P4ll35l . 

Discrete Hardy space h 1 . Take a smooth, compactly supported cut-off function 
v supported in the centered unit cube Cq, with v = 1 near the origin. Then define 
a "discrete Riesz kernel" sequence z = {zfc} fcgZ d G I 2 by specifying its Fourier 
series: 

C(£) = £ z k e~ 2 ^ k := -i-^-MO, £ G Co, (5) 

where £k denotes the dot product (recall £ is a row and A; is a column vector) 
and where for later convenience we use — £ rather than +£ in the exponent of the 
Fourier series. 

The sequence z is vector-valued (since z k G C d ), and belongs to £ 2 because 
(— z£6 _1 /|£6 _1 |)i/(£) is bounded and hence belongs to L 2 (Cq). Thus the series for 
the periodic function ( converges unconditionally in L 2 {Cq). When b = I, observe 
from ((5]> that £ is simply a cut-off version of the Fourier transform of the Riesz 
kernel. 

Naturally s * z G £ 2 whenever s G I 1 . We define a "discrete Hardy space" by 
requiring that s * z belong to the smaller space l x : 

h 1 = {s G i 1 : s*z G l 1 } , 

with a norm 

ll s lki : = INI^ 1 + II s * zh 1 

that makes h 1 a Banach space. 

Appendix IBIinvestigates some properties of h l , including its independence from 
the cut-off function v, and its relation to the atomic sequence space B\Z d ) stud- 
ied by several authors. The appendix also points out that sequences in h l have 
vanishing mean: J] fegZ d Sfe = 0. 



with the norm := ||/||i + \\Rf\\ v 
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The mixed-norm Banach space we will need is 

£ 1 (h 1 ) = {ce£ 1 {£ 1 ):c*ze£ 1 {£ 1 )} 
= {ce£ l {£ 1 ) : IMI^i) < oo}, 
where the convolution is taken with respect to the fc-index and the norm is 

INI^/i 1 ) := llcll^ 1 ^ 1 ) + H c * z \\^{^) 

j>0 j>0 k&L d 

with the notation cj = {cj t k}keZ d - 

Properties of the synthesis operator. We will prove boundedness and surjectivity 
of the synthesis operator. Note 

^j,k( x ) = I det a,j\ip(a,jX — bk) 

in the next two theorems, because we implicitly take p = 1 when working with 
ip £ L l . 

First we show boundedness. 

Theorem 7 (Synthesis into H 1 ). Assume ip € L 1 and 

sup \\ip - ip{- -y)\\m < °o. (6) 

lvl<i 

Then S : ^(h 1 ) — > H 1 is bounded. More precisely, if c £ £ (/i ) the series 
Sc = J2j>o J2k& d c i,k4>j,k converges in H 1 in the sense that 

the sum over k converges absolutely in L 1 to a function belonging to H 1 
and the sum over j converges absolutely in H , 

and furthermore \\Sc\\jji < C||c||^i(/ji) for some constant C = C(ip,b). 

See Section [9Tl for the proof. To understand why we so carefully describe the 
convergence of Sc in H l , in (0, just recall that ipj^ ^ H 1 when f Rd ip dx ^ 0. 

Many functions ip satisfy the finite supremum assumption in ©, for example if 
ip € L p for some p > 1 and ip has compact support, or if ip is a Schwartz function; 
cf. JU §3.3-3.4]. Incidentally, the supremum in (O can equivalently be taken over 
any other ball of y- values, by (H Lemma 6]. 

Theorem [7j was proved for synthesizers ip € L 2 having compact support by S. 
Boza and M. J. Carro in Proposition 3.11], [3] Theorem 3.1]. Their methods 
are very different from ours, involving a maximal characterization of H p (Ij d ),0 < 
p < 1. Theorem [7j was also proved earlier in one dimension for synthesizers 
ip € having compact support and locally integrable Hilbert transform Rip € 

Lj oc (M.), by Q. Sun 071 Theorem 13]. Note that Sun's assumptions on ip imply 
the condition ©. Interestingly, Sun proved a converse theorem for compactly 
supported ip, saying that bounded synthesis implies Rip must be locally integrable. 

Now we show surjectivity of synthesis from i 1 (h 1 ) to H , a result that seems to 
be qualitatively new. 
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Theorem 8 (Synthesis onto H 1 ). Assume ip G L 1 with J Rd ip dx = 1 and 

H - - y)\\m ^ asy^O. (8) 

Then S : £ (h 1 ) — * H 1 is open, and surjective. Indeed if f G H 1 and e > 
then a sequence c G ^(h 1 ) exists such that Sc = f (with convergence as in 
and 

||c||£i(hi) < C\\f\\ H i +e, 
for some constant C = C(b). 

See Section [9]4] for the proof. The constant C can be evaluated, if desired. 

Assumption ([U) holds, for example, if ip G -L p for some p > 1 and ^ has 
compact support, or if ip is a Schwartz function (cf. (U §3.3]). 

By combining the last two theorems, we obtain a norm for H 1 in terms just of 
coefficients in affine expansions. 

Corollary 9 (Synthesis norm for H 1 ). Assume ip G L 1 with L d ipdx = 1, and 
suppose ([8]) holds. Then for all f G i? 1 , 

w inf { Hcll^i^i) : c G £ (h ) f = Sc as in (0} . 

The proof is easy: assumption ([U) implies © by [6, §3.4], and so Theorems [7] 
and[8]both apply here, giving the corollary. 

The analysis operator also provides a coefficient norm for H 1 : 

Corollary 10 (Analysis norm for H 1 ). Take an analyzer <p G L 1 with f Rd (pdx = 

1, and suppose (p is a Schwartz function with (p supported in C^b^ 1 . 
Then for all f G H\ 

WfWsi « \\Tf\\^ m =sup||{(/,^, Jfe )}|| h i |detfe|. 

The corollary says the analysis operator T is linear, bounded and injective from 
H 1 onto its range in ^(h 1 ). See Section 1931 for the proof. 

5. Sobolev space results 

Write W m,p = W m,p (W i ) for the class of Sobolev functions with m derivatives 
in LP, normed by 

11/11*-* = £ Wf\\ v . 

\p\<m 

Here p = (pi, . . . , p^) is a multiindex of order \p\ = pi + • • • + Pd- 

We continue to assume the dilation matrices are isotropic and expanding, with 

a,j = djl, like in the previous section. 

We first construct a sequence space on which the synthesis operator will act, and 

then construct a class of synthesizers, before stating our Sobolev synthesis results. 

Discrete Sobolev space. Define difference operators on sequences s = {sk\k& d 

by 

A t s = {s k - s k - et } k& d, 
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for t = 1, . . . , d, where e t is the unit vector in the t-th coordinate direction. Define 
higher difference operators by 

A p s = Af • • • A p /s. 
Then we can define a discrete Sobolev space by 

w m ' p = {s G £ p : A p s G £ p for each multiindex p of order \p\ < m} , 
with norm 

\p\<m 

This space is just £ p with a new norm, but we proceed to weight the norm by 
appropriate powers of the dilation sequence, as follows. 

Recall sequences are multiplied term-by-term, meaning (cc)j j k = Cj^Cj^k and 
so on. In particular we can define powers a r = {a^}j>o of the dilation sequence, 
whenever r is a nonnegative integer. With these conventions, we define a dilation- 
weighted discrete Sobolev space by 

£V m,P , «) = {c G £ l (£ p ) : a^A p c G £ l {f) for each multiindex p of order \p\ < 
with norm 

||c||^i( w m,P, a ) := ||a' p 'A p c||^i(£p) 
|p| < m 

i>o |p|<m \fcez d 

(The difference operators here should be understood as acting on the fc-index of the 
sequence Cj One can check £ 1 (w m,p , a) is a Banach space. 

77je class of synthesizers. Our synthesis results will hold when ip has the special 
convolution form 

m factors 

, * s 

ip = /3*---*p*rj, (9) 

where 

p=\bc\- 1 t bC 

is the normalized indicator function of the period box in the lattice bZ d . (One can 
show that the convolution form (© is equivalent in one dimension to a Strang-Fix 
condition on ip, but it is definitely stronger in higher dimensions, as explained in 
(7J Notes on Theorem 1].) The point of this convolution form © is that derivatives 
of tjj turn into differences, in formula (|23T ) later on, and these differences transfer 
to the coefficient sequence in formula ( |27l) . 

Now we can state the boundedness of the synthesis operator on Sobolev space. 
Our statement involves the matrix rescaling operator 

(M b f)(x) = \detb\f(bx). 

Choosing b = I gives the simplest results in what follows, of course. 
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Theorem 11 (Synthesis into W m ' p ). Assume 1 < p < oo and rj G L p with P\rj\ G 
Lf oc . Le£ m G N awcf define ip by ©. 

ip G W m ' p , and S : £ 1 (w m ' p ,a) — » VF m ' p jj bounded. More precisely, if 
c G £ 1 (w m ' p , a) f/ierc senes 5c = X^>o ^kez d c j,k^j,h converges in W m ' p in 
the sense that 

the sum over k converges pointwise absolutely a.e. to a function in W m ' p _ 
and the sum over j converges absolutely in W m,p , 

and furthermore 

\\M b Sc\\ W m, P < |det6|" 1/p ||P|?7||| L p( &C )||c||£i( lu m, P>a ). 

See Section [10. II for the proof. Prior results on bounded synthesis include ll28l 
37 ], for compactly supported ip. Here tp need not be compactly supported. 

Regarding the appearance of the rescaling operator Mb in the theorem, note the 
analogous L p inequality ® in Theorem [JJ can be put in the same form as above, 
namely 

\\M b Sc\\ p < \ detb\- 1 ^PM\ IJ , {bC) \\c\\ eHeP) , 
by applying Theorem [T] and then making a change of variable on the lefthand side. 
Now we obtain surjectivity of S. 

Theorem 12 (Synthesis onto W m,p ). Assume 1 < p < oo and rj G LP with 
P\rj\ G L p oc and f Rd rj dx = 1. Let m G N and define ip by ©. 

Then S : i l {w m < p ,a) — > W m ' p is open, and surjective. Indeed if f G W m < p 
and e > then a sequence c G l 1 (w m ' p , a) exists such that Sc = M^ 1 f with 
convergence as in ( fTOl) , and such that 

\\ c \\e 1 (w m -P,a) - ll/llw™* + e - 

The theorem is proved in Section 110.41 Note the conclusion of the LP result 
Theorem [2] can be rephrased to look like Theorem [T2l by applying it to M h ~ l f 
instead of to /. 

There seem to be no direct predecessors in the literature for the surjectivity in 
Theorem [12] Indirect predecessors include the density results from Strang-Fix 
theory, which we discuss in Q §3.5] and after Theorem l24l below. To learn about 
wavelet and affine frame expansions in Sobolev space, using all scales j G Z, 
consult ll29l l32ll and the recent work in iflOl . 

Norm equivalence now follows from Theorems ITTI and [T2l 

Corollary 13 (Synthesis norm for W m,p ). Assume 1 < p < oo and rj G L p with 
P\rj\ G L p oc and j Rd rj dx = 1. Let m G N and define ip by ©. 
Then 

\\f\\ W m, P Kt inf {\\c\\ e i {wm ,p ta) : M b ~ l f = Sc as in £[0]>} 
for all f G W m ' p . Explicitly, 

|det6| 1 /f||P|, 7 |||2 p 1 (6C) ||/|| wm>P <inf{||c||^i ( ^, Pja) :M b - 1 f = Sc as in m} 

< \\f\\w m <p- 
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The coefficient norm can equal the standard Sobolev norm: 

Corollary 14 (Synthesis norm equality). Assume 1 < p < oo and that rj S LP is 
nonnegative. When p = 1 assume J^ d r/dx = 1, and when 1 < p < oo assume 
Prj = 1 ( which implies J Rd r/dx = 1). Let m £ N and define by ([9]). 
Then for all f (£W m > p , 

\\f\\w m -p = inf { Hcllp^m.p^) : M^ x f = Sc as in (PTOl)} . 

The analysis operator gives a norm to Sobolev space also: 

Corollary 15 (Analysis norm for W m ' p ). Assume 1 < p < oo and m G N. Ta&e 
an analyzer (j) £ L q with P\(f>\ € LP° and L d cf>dx = 1, and assume b = L. 
Then for all f eW m ' p , 

\\f\\w™<p ~ \\Tf\Uoor w m, Pj(x) = sup V lajl^HATj/H^ 

i>o , ~~ 

|p|<m 

= sup V |a,|H||{A^/,^ fc )}|| £P . 
i>o , *r? 

\p\<m 

Thus the analysis operator is linear, bounded and injective from W m,p onto its 
range. Section [103] has the proof. 

6. Connection to phi-transforms and Triebel-Lizorkin spaces 

The decomposition or representation of function spaces by means of an affine 
system generated by a single function (or collection of functions) is a well-established 
technique in harmonic analysis. Important examples include wavelet expansions 
and the phi-transforms. For wavelet theory we refer to (H [13] [29] [32]] and the ref- 
erences therein. The connection between phi-transform theory and our results is 
briefly explained below. 

Triebel-Lizorkin spaces include LP ', H 1 and W m ' p , 1 < p < oo, but neither L 
nor VP" 1 ' 1 . Phi transform theory implies the following results about these spaces; 
see HI [22] [23] EH for a complete account. 

(i) That the Triebel-Lizorkin space norm of an element / is equivalent to the in- 
fimum of the corresponding sequence space norm of c, where the infimum is taken 
over all representations f = Sc and the synthesizing tp satisfies some moment 
condition (in particular J Rd dx = 0), and a decay condition and a "Tauberian" 
condition. Note that the dilations are assumed isotropic, the translation matrix b 
must be "sufficiently small", and that the sequence space norm in this theory also 
involves an integration with respect to the continuous variable on W 1 . 

(ii) That there exists an analyzing function <p (or a collection of functions 0^ ,fe )) 
satisfying similar conditions to tp such that the Triebel-Lizorkin space norm of / 
is equivalent to the sequence space norm of Tf = {(/, <pj &)}, and / is represented 

by/ >;,,,*//;,/•. 

The large scales j < must be included in the synthesis (i) and the analysis (ii). 

Versions of these results have also been proved for affine systems generated 
by interesting families of functions (instead of one function ip); see |[38l and the 
references there for "Gausslet" and "Quarkonial" analysis. 
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Our main results in Sections[3]-[5]can be viewed as analogues of the above results 
when both the analyzer <p and the synthesizer ip have non-zero integral (that is, no 
vanishing moments). Further, our sequence space norm is simpler, involving only 
discrete sums and no integration with respect to the continuous variable. And the 
norm equivalence statements in our Corollaries [6l \\0\ and [15] take an especially 
simple form. These corollaries are in the spirit of norm equivalences in frame 
theory, but in our case the frame decomposition / = J2j k SjTjf cannot hold, 
because if / i— > Ylj k SjTjf is bounded on L 2 , then either <fi or ip must have a 
vanishing moment by |[25l Theorem B]. Nevertheless, our corollaries give rise to 
Banach frames in the sense of K. Grochenig ll26l . Details are given in Appendix 



7. Open problems 

Throughout this paper we have assumed the synthesizer ip has nonvanishing 
integral. When it has vanishing integral, J Rd ipdx = 0, we do not have a compre- 
hensive understanding of conditions under which one can synthesize surjectively 
onto L p or Hardy or Sobolev space. Of course there are substantial classes of syn- 
thesizers ip for which surjectivity and even injectivity is known, for example the 
class of wavelets lfl"3l 12911321 (provided one includes large scales j < in the affine 
systems). But wavelets seem rather restricted objects to us. In accordance with the 
goal expressed in the Introduction, we conjecture that some much more general 
surjectivity result should hold when ip has vanishing integral. Recent L 2 work of 
Gilbert et al. |[25l Theorem G] is a step in the right direction, for it assumes only 
that ip(0) = and ip has some cancellation properties. Unfortunately the result 
suffers from oversampling of both dilations and translations, and ways to remove 
that oversampling remain a mystery. 

Problems with vanishing integrals can be more challenging than they appear. 
For instance, it is an open problem of Y. Meyer 021 p. 137] to determine whether 
the affine system {ip{2^x — k) : j,k G Z} spans L P (M) for each 1 < p < oo, when 
ip(x) = (1 — x 2 )e~ x I 2 is the Mexican hat function (the second derivative of the 
Gaussian). This is known to be true when p = 2, since the system forms a frame, 
but it remains open for all other p-values. It is apparently also open to determine 
whether the system spans H 1 (R). To express this in terms of the synthesis operator, 
first notice S$ is bounded from ^(l 1 ) to H 1 because the Mexican hat ip belongs 
to H 1 , and then ask: is : £ (£ l ) — > surjective? Surjectivity would give 

an atomic decomposition of H 1 in terms of the Mexican hat affine system. 

There are two partial results in the literature dealing with the Mexican hat prob- 
lem. In ll24l the authors proved that there exist (sufficiently small) r > 1 and 
s > such that the affine system {ip(r^x — sk) : j,k G Z} spans H P (R), 
1 < p < oo, while in (H the authors proved the same result for the affine system 
{ip{2^x — bk) : j,k G Z}, where b > is sufficiently small and 1/2 < p < oo. 
Since H' p = LP for 1 < p < oo, these results show that the spanning property 
holds for the Mexican hat function provided we accept some degree of oversam- 
pling. 
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We conclude by pointing out a gap in understanding of our Hardy space synthe- 
sis results. We do not know whether h 1 is the "natural" domain for the synthesis 
operator in Theorem [7] Can one prove it is natural, in the sense that 

J, s k Tp jtk G H 1 => sell 1 

k& d 

whenever tp G L 1 satisfies hypothesis ©? 

8. W proofs 

8.1. Proof of Theorem [J— synthesis l x ($P) -» LP. 

First assume 1 < p < oo and ip G LP with P\ip\ G LP, . We will synthesize at a 
fixed scale j > 0, by taking s£f and defining 

/ = Sjs = ^2 s k ipj,k- 

k£Z d 

The task is to show / G LP with ||/|| p < ||s||^p||P|V'|||lp(6C)/I det 6|. Then Theo- 
rem [Qfoilows easily, by summing over the dilation scales. 
We have 

{ V 

I s k ipj,k{x)\ p < |sfc||V(aj-ar — 6fe)| | detail 



— I^H^C^i 37 — ^)| IVK .?^ — bk)\ I | detail 

by Holder's inequality on the sum. Integrating with respect to x yields that 

II s ^M\ P p < E \ S *\ P f ~ bk)\(m(x)) p - 1 dx/\ det b\*- 1 

fcgZ d k&L d Rd 

= | Sfc | p ||P|^|||i P(6C) /|det6r (11) 

by changing x ^ x + bk and then periodizing the integral. We conclude that the 
sum over k in (Q]) converges pointwise absolutely a.e. to an L p function, and that 



E s k^j,k\\p ^ ll s IM|PMIlLP(fcC)/|det6|. 



It remains to prove the theorem when p = oo. So assume p = oo and ip G L°° 
with P\ip\ G L^ c . Then P|*0| is bounded, since it is locally bounded and periodic. 
If c G e l (£°°) then 

I E E ^fc^fc^)! < E( sup l c i. fe D E IVK ^ - bk )\ 

j>o kez d j>o kGZ,d k ez d 

< ||c||^i,oo llP^llloo/l det 6| < oo 

for almost every x, from which the theorem follows. 
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Is £ 1 (£ p ) the correct domain for synthesis? 

We have proved S is bounded from £ 1 {£ P ) into L p . Could S be bounded on an 
even larger domain? The natural candidate would be £ r (£ p ) with r > l,p' > p, 
but we will show by example that 5 need not be bounded on this domain unless 
r = 1 and p' = p. 

Work in one dimension with 6=1 and dyadic dilations aj = 2 J , and choose 
to be supported in the unit interval [0, 1). Then for any sequence s G £ p ' we have 
Sjs(x) = J2 keZ Sk2^ p ip(2 : > x — k), which has L p norm ||<Sjs|| p = H^llpIN^- 
Hence if SjS G L p then s € £ p , so that for S to be bounded on £ r (£ p ') it is 
necessary that p' = p. 

Further, if t = {tj}j>o is an y nonnegative sequence in £ r then the sequence 

it i 2~i/ p , k = 0,l,...,Z> -1, 

Cjk= \ 

I 0, otherwise, 

belongs to £ r (£ p ), and if tp = lp i) is the indicator function of the unit interval 
then Sc = \\t\\£i l[o,i)- Thus for Sc to belong to L p it is necessary that t G £ 1 . 
Hence r = 1, as claimed. 

So we cannot expect to enlarge the domain £ 1 {£ P ), in general, when f R tp dx ^ 
0. But there is a loophole relevant to wavelets, because if f^tpdx = then the 
t/jj ; k might exhibit some cancellation between different j-scales that invalidates our 
"t" example above and allows us to take r > 1. For example, if p = 2 and ^ is a 
wavelet (meaning the functions ^ £ are orthonormal and complete in L 2 ) then S is 
not only bounded but is an isometry from £ 2 (Z x Z) to L 2 . Thus the natural domain 
in the wavelet case has r = 2. Recall integrable wavelets satisfy f R ip dx = by 
lUll p. 348]. 

Continuing in the special case p = 2, the ipj^ form a frame (which is more 
general than an orthonormal basis) if and only if the synthesis operator is bounded 
and surjective from £ 2 (Z x Z) to I? . This is a special case of Christensen's Hilbert 
space result [9, Theorem 5.5.1]. Thus the natural domain for frame synthesis has 
r = 2. This agrees with our earlier remarks, because integrable frame generators 
are known to satisfy J R ip dx = 0. 

8.2. Analysis L p -> £ p . 

The proof of Theorem |2] relies on the following estimate for analyzers <p with 
bounded periodization. Recall the analysis operator Tj at scale j was defined in 
Section 12 

Proposition 16 (Analysis into £ p ). Assume 1 < p < oo and <j) G L q with P\<p\ G 
L°°. Then for each j , 

Tj :L P ^£ P with norm \\Tj\\ < | det H^. 

The proposition is known from Aldroubi, Sun and Tang [1 , formula (2.2)]. 
The hypothesis that \<fi\ have bounded periodization means <fi is bounded and its 
integer translates "do not overlap too often". 
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Proof of Proposition [i6l Let / G LP. When 1 < p < oo we have 

WTjfWtP 




by Holder's inequality on the inner product 



= 1 det b\^ (J \mppMw) dy) 1P U\\l /q 

< Idetfol^llP^HUII/llp, 

using that ||^||i = f K Y,k& d \^( x ~ bk )\ dx ^ ll-P^I IU- 

When p = 00, the proof is straightforward. □ 

Aside. For the special case p = 2, Proposition [T6l is known ||9l Theorem 7.2.3] 
for all 4> G L 2 with P{\4>\ 2 ) G L°°, and one can show this condition is weaker 
than P\4>\ € L°°. That is, analysis is bounded from L 2 to £ 2 with norm ||Tj|| < 

||-P(|0| 2 )||oo 2 > where the periodization is with respect to the lattice Z d 6 _1 . This 
is another way of saying the translates of (f> form a Bessel sequence, or satisfy an 
upper frame bound. 

8.3. Scale-averaged approximation in LP. 

The following approximation result will be used in proving Theorem [2] (surjec- 
tivity of the synthesis operator). The result is interesting in its own right too, due 
to its explicit nature: we simply analyze with <fi, then synthesize with ip, and then 
average over all dilation scales to recover /. 

Say that the dilations expand exponentially if 

ll a i a 7+iH — ^ for all j > 0, 
for some < 5 < 1. In one dimension, this means \a,j \ is a lacunary sequence. 

Theorem 17 ([5, Theorem 1 and Lemma 2]). Assume 1 < p < 00 and tp € 

LP,P\ip\ G Lf oc ,0 G L",P\4>\ G / G LP. Assume f Rd 4>dx = I, and write 
1 = J Rd ip dx. 

(a) [Constant periodization] If Pip = 7 a.e. then 

SjTjf — > 7/ in LP as j — > 00. 

(b) [Scale-averaged approximation] If the dilations aj expand exponentially, then 

1 

— SjTjf — > 7/ in LP as J — > 00. 
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The scale averaged approximation in part (b) can be written in full as 



±£|det&|[ 
\ 



Y (/' 0j,k)ipj,k ) -> 7/ in LP > as J ^ 00. 

fcez d 



Aside. Part (a) of the theorem has a long history, summarized in [5J §3]. 

8.4. Proof of Theorem|2] — synthesis onto LP. 

We can assume the dilations expand exponentially, as follows. For each j we 
have || ajaj +r \\ < \\ctj\\ \\aj +r \\ — > as r — > 00 because the dilations are expanding. 

Thus llajaj+rll < 1/2 provided we choose r sufficiently large. By iterating this 
argument we arrive at a subsequence of dilations that expands exponentially (with 
5 = 1/2). It is enough to use only this subsequence of dilations, when proving 
openness and surjectivity of the synthesis operator. 

Consider / G L p and J G N, and define a sequence cj = {c j;j,k}j>o kez d by 

cj. jk = hdetb\f {f ^ k) «*j = l,...,J, (12) 
J 10 otherwise, 

where = |6C| _1 l(,c is a normalized indicator function. Note P\<f>\ = 1. Then 
cj G t x (l v ) because applying Proposition [TBIfor each j = 1, . . . , J gives that 

1 J 

IM| w < - Y, I det 6| 1A 1/|| P = I det 
J i=i 

And clearly 

J 



Scj = j^2s j T j f 



— > / in L p as J — > 00, 

by Theorem fTTl b). (Here we use that the a,j expand exponentially.) Thus the open 
mapping theorem in Appendix IA1 says S : — s> L p is open and surjective, 
and that for each / G L p and e > there exists c G i 1 ^) with Sc = f and 
||c|| w < Idetb^WfWp + e. 

8.5. Proof of CorollaryU] — Synthesis norm for L p (f2). 

The "<" direction of the norm equivalence follows immediately from Corol- 
lary 13 since we are restricting the collection of sequences c that can be used to 
represent /. 

For the ">" direction, we will modify the proof of Theorem|2] Like in that theo- 
rem we can assume the dilations expand exponentially, by passing to a subsequence 
of j-values. 

Suppose / G L P (Q) is supported on a compact subset of O. We claim there 
exists jo > such that 

spt#j,fc) C whenever j > j and (/, (f> jik ) ^ 0, (13) 
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where <fi = |6C| _1 l&e as previously. The existence of j should be clear intuitively, 
since the support of / lies at some positive distance from the boundary of Q while 
tp and 4> have compact support and WaJ 1 1| — > 0. We leave the detailed proof of (fT3T > 
to the reader. 

Next, consider J G N and define 



c J;j,k 




for j =j + l,...,io + J, 
otherwise, 



-I JU I " 

cj\\p m <- y, i det6 i 1/, n/iip = i det6 i 1/ *n/iip- ( 14 ) 



so that cj G £ x (£ p ) by Proposition [T6l with 

jo+J 

J . 

3=30+1 

We know Scj — > / in L p (R rf ) as J — > oo, by Theorem [TTlb) with an index shift 
on the dilations. 
Define 

£ = {c G i l (i p ) : c is adapted to fl and ip}. 

Clearly £ is a closed subspace of I 1 (£ p ), and hence is a Banach space. Note cj £ £ 
by {□]), so that 5c = on ft c . 

We have verified the hypotheses of the open mapping theorem in Appendix |A] 
for S : C ->■ L p (f2), with the constant ^4 = | det 6| 1/<? by £[4]>. Admittedly we 
have verified the hypotheses only for the dense class of / having compact support 
in Q, but a dense class suffices, by the comment at the end of Appendix |Al The 
open mapping theorem tells us S : C — > L P (Q) is surjective, and that for each / G 
L p (fl) and e > there exists c G C with Sc = f and ||c||£i^ P ) < | det &| 1 ^ 9 ||/||p + 
e. This proves the ">" direction of the corollary. 

8.6. Proof of Corollary [6] — analysis norm for L p . 

By boundedness of the analysis operator in Proposition [T6j 

sup < Idetftl^HPHIIooll/llp. 

3 

To prove the other direction of the norm equivalence, choose a function ip that 
satisfies the hypotheses of Theorem [TTl a) with 7=1. Then by that theorem, 
SjTjf — > / in LP as j —> 00. Therefore it follows from Theorem [T] (bounded 
synthesis) that 

p < sup \\SjTjf\\p < |det6| _1 ||P|^|||LP(6c) SU P \\ T jfh p i 
3 3 



which proves the corollary. 



9. Hardy space proofs 



In this section we assume the dilations are isotropic and expanding, so that a 



.1 



ctjl for some nonzero real numbers otj with \aj\ — > 00. 
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9.1. Proof of Theorem|7]— synthesis P-Qi 1 ) — ► H 1 . 

We begin by showing that the Riesz transform "almost" commutes with affine 
synthesis, which will be the key step in proving boundedness of the synthesis op- 
erator, in Theorem [7] 

Lemma 18 (Riesz transforms commute with synthesis). Suppose ip E L 1 . Take v 
to be the smooth, compactly supported cut-off function used to define z in formula 
©. Let p, be the Schwartz, function with = v{£,b), and let X be a Schwartz 
function with A supported in C^b^ 1 and with A(0) = 1. 
If s E h 1 then 

Z * ( E Sk(ip * A * fi)(x — bk)) = (z * s)k(ip * A)(rr — bk) 
k£Z d kez d 

EL 1 . 

We can rephrase the lemma (after rescaling x ajX and using the dilation 
invariance of the Riesz transform) as saying 

RSj^u^s = sign(a j)Sj^\Rs, s E h , 

where on the lefthand side R denotes the continuous Riesz transform (convolution 
with the Riesz kernel Z) and on the righthand side R denotes the discrete Riesz 
transform (convolution with the discrete Riesz kernel z). Thus we see Lemma [T8l 
is a discrete analogue of the formula R(ip * f) = ip * (Rf), once we remember that 
affine synthesis is a discrete analogue of convolution with a synthesizer. 

Proof of Lemma U8\ Define f(x) = J2kez d s h{i ) * ^ * ~~ This sum 
converges absolutely in L , because 

ll/lli < E \ s k\U* A*mO = ||fl||^||^*A*A*lli- 

Our task is to show / E H , with its Riesz transform being as stated in the lemma. 
Consider the periodic functions 

*(o = E ^ e ~ 2 ^ fc < ceo = E ^ e "* rfe *» 

which are well defined since sGl 1 and z£l 2 . We have 

= -i||^(^)^)A(0?(0 

= C(#MW(WO d5) 

by definition of the Riesz kernel sequence z in ([5]>, using here that A(£) = when 

Of course C(£) a (£) = Sfc6Z d ( z * s )k^~ 2m ^ k m L 2 (Co), by computing Fourier 
coefficients of the two sides in L 2 (Cq) and using z E £ 2 ,s E ^,2 * s E ^ 2 . 
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Therefore (fT5l) says 

— «-j7j-/(£) = Fourier transform of ^ (z * s)fc(*0 * A) (a; — &&), 

where we note that z * s G -f 1 by the hypothesis s G ft 1 . Thus (Z * /)(x) = 
Sfeez^l- 2 * s )fc(V ; * A)(x — 6/c) G L 1 . This proves the lemma. □ 

Proof of Theorem\7\ Fix j > and take s G ft, 1 . Define 

/ = ^ Skipj^. 

k£Z d 

We will show that / G H 1 with < C\\s\\ h i, where C = C(ip,b) is inde- 

pendent of j. Then Theorem [7] follows by summing over j. 
It is enough to show that the function 

g(x) = (M-/f)(x) = £ ^{x - bk) 
k&z d 

belongs to H 1 with ||<?||#i < C||s||^i, because = IIsIIh 1 (using that 

aj = djl is isotropic). Obviously \\g\\i < \\s\\ei and so our task is to show 

\\Rg\\i < C\\s\\ h i. 

To understand Rg, take u, \i and A as in Lemma [HI and decompose 
9 = 9i + 92 where 

9i{x) = ^2 Sk ^ * A * ^ x ~ bk ^ 

kez d 

92{x) = ^2 Sk ^ - ^ * A * fi)(x - bk). 

k&L d 

Lemma [18] implies 

ll-Rffilli < \\ z * s\\e-\\ip * A||i < * A||i. 

(Thus we see Lemma [18] is used to push the Riesz transform onto the coefficient 
sequence s, which belongs to ft 1 , rather than onto the synthesizer i/j * X* fi, which 
does not belong to H 1 .) 

The sum defining gi converges absolutely in H , because s G ft 1 C i 1 by 
assumption and ip — ij) * \ * n £ H l by Lemma [l9]below (which is where we use 
the hypothesis (O). Hence 

\\R92\h < INI^H^ - i> * A * A*) ||i. 
These bounds prove Theorem [7] 

We must still prove Lemma[[9] needed to treat g 2 in the "smoothing" step of the 
proof above. We use the translation operator Tyifi = ip(- — y). 
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Lemma 19. Suppose ifi G L 1 and \\R(ip — T y ip)\\i < C < oo for all \y\ < 1, and 
A G L 1 w;Y/j |y|A(y) G L 1 OM<i J Rd A(y) = 1. Then 

\\R{il)-il>*\)h<C [ (\y\ + l)\X(y)\dy. 

JR d 

The assumption that \\R(ip — T y ip)\\i < C for \y\ < 1 is a restatement of hy- 
pothesis ([6]). 

Proof of Lemma[[9\ The first step is to show the Hardy norm of a difference grows 
at most linearly with the difference step, that is 

\\R^-Tyi>)\W<C{\y\ + l) forallyGM d (16) 

For this, let m be the integer satisfying \y\ < m < \y\ + 1. After writing y as a 
sum of m vectors each having norm less than 1, we can prove (fTBT) with Cm on 
the righthand side by telescoping the differences and using the triangle inequality, 
noting the Riesz transform is translation invariant. 
Now observe the function 

ip(x) - (ip* X)(x) = / (tp(x) - T y ip(x))X(y) dy 
belongs to H 1 and has Riesz transform 

R(i> - tP * X)(x) = [ R(il>-Tyfl))(x)\(y)dy, 

by El Lemma 10]. That is, one can take the Riesz transform through the integral. 
Thus 

|W-V*A)[|i< / WR^-Ty^UX^dy 

JR d 

< [ C(\y\ + l)\X(y)\dy 

JR d 

by ( fT6l ), as desired. □ 

9.2. Analysis H 1 —> h 1 . 

For the proof of Theorem [8] we want boundedness of the analysis operator Tj 
from H 1 to h 1 , which we state in Proposition |2"T1 below. First we show: 

Lemma 20 (Riesz transforms commute with analysis). Assume the analyzer cj) is 
a Schwartz function with (f) supported in Cob^ 1 . Fix j > 0. Choose v as in the 
definition of h 1 in Section [?] and let fij be the Schwartz function with /U^(£) = 

IffeH 1 then 

RTjf = sign(a j )T i ii(/x j * /), 
where on the lefthand side R denotes the discrete Riesz. transform ( convolution with 
the discrete Riesz. kernel z) and on the righthand side R denotes the continuous 
Riesz. transform (convolution with the Riesz kernel Z). 

Remember ctj denotes the isotropic dilation factor in a-, = cijl. 
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Proof of Lemma \20\ Observe 4>j^(x) = 4>(ajX — bk) in what follows, because we 
implicitly assume p = 1, q = oo, wherever we deal with the Hardy space. 
The k-th term of the sequence RTjf is 

(z * Tjf) k = | det b\ ^2 ^(/' 

= | det b| £ f zee-^fteajMSe-^dt, 

eeZd JC b 1 

by Plancherel and the compact support of 4>. Substituting in the definition of ( 
from ([5]) and then using the definition of fj,j, we find 

(* * Tjf) k 

= \detb\[ ^■v(&)f(Za j )4(0e-W<dZ 
Jr* Isl 

= | det b\ / siga(aj)(fj,j * RfT(€aj)<j)(€)e- 2 ™& k dt; since Oj = ay/ 

= | det 6|sign(a i ) *f),<f>j,k) 
by Plancherel, and this is the fc-th term of sign(a>j)TjR([ij * /), as desired. □ 

Proposition 21 (Analysis into h l ). Take cj> and v as in Lemma \20\ Then for each 
3, 

Tj : H 1 — » h 1 with norm \\Tj\\ < \\P\<fi\ ||oo ll^lli- 
Proof of Proposition \2l\ If / G H 1 then 

\\TjfU < ||P|0|||oc||/||l 

by Proposition [TBI with p = 1, while 

11**7}/^ = \\RTjf\\ el < ||P|0[||oo||i2(/*i*/)||i 

< l|P|0Hloo|KI|l||P/||l 

by combining Lemmal20land Proposition [16] Add these two estimates and observe 
ll/^jlli = ll^lli — 1^(0)1 = 1. by definition of fij in Lemma l20l □ 

Aside. Our compact support assumption on the Fourier transform of the analyzer 4>, 
in Lemma|20]and Proposition [21] seems rather strong. Perhaps it can be weakened. 
But notice it ensures that ^(Ib^ 1 ) = for all row vectors £ G Z d \ {0}, which 
implies P(f> = const. This constant periodization condition is necessary, as follows. 
If Tj maps H 1 into h , then for all / G H 1 we have 

/ f(x)P<f>( aj x)dx = V (T 3 f) k = 

by the zero-mean property of Tjf G h 1 (see Appendix [Bj. Taking / G H l to 
approach a difference of delta functions implies Pcj) is constant. 
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9.3. Scale-averaged approximation in H 1 . 

In this section we prove scale-averaged convergence in H , which we need for 
the proof of Theorem [U 

Theorem 22. Assume ip G L 1 with L d ip dx = 1 and 

\\ip - - y)\\m ^ asy^O, (17) 

which is hypothesis ([8]). Assume <f> is a Schwartz function with <f> supported in Co^ 1 

and f Rd (j) dx = 1. Lef / G JET 1 . 

(a) [Constant periodization] If Pip = 1 a.e. then 

SjTjf ~^ f in H l as j — > oo. 
[Scale-averaged approximation] dilations a>j expand exponentially, then 

1 J 

— SjTjf — > f in H 1 as J — > oo. 
i=i 

The theorem was proved in our earlier paper [6, Theorem 1] for analyzers (f> with 
compact support and P<j) = const., by comparing SjTjf with an approximate iden- 
tity formula. The new proof below is more conceptually satisfying, as it is based 
on commuting the Riesz transform through the analysis and synthesis operators. 

Proof of Theorem 122] Convergence in L , for parts (a) and (b), follows immedi- 
ately from our TP result Theorem [T7J with p = 1. 
To prove convergence in H , we want to show 

RSjTjf -> Rf in L 1 , for part (a), and (18) 

1 J 

- V i25jTj/ -> i?/ in L 1 , for part (b). (19) 

To begin, suppose = z^(£&) and A are as in Lemma [T8l and decompose 

SjTjf — j&Tjf — Sj,ij}—'ip*A*tJ,Tjf -f- Sj^xXx^Tjf 
= Aj + i?j , say, 

where we write Sj^ and so on to emphasize the synthesizer being used, in each 
part of the formula. We have 

RBj = sign.(atj)Sj^*\RTjf by Lemma [T8l 

= Sj^*\TjR(p,j * f) by Lemma|20] 

= Sj^*\Tj([ij * (Rf) — Rf) + Sj^xTjRf 

= Cj + Dj , say. 

We estimate Cj by 

||Cj||i — (const.) * (Rf) — Rf\\i by Theorem [T1 and Proposition [T6l 

— > as j — > oo, 
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since Hj(x) = | det b~ 1 aj\u(b~ 1 ajx) is an approximation to the identity (recalling 
aj = ctjl with \aj\ — » oo). 

Theorem [T7] implies Dj — > Rf in L 1 in part (a), and implies in part (b) that 

Thus to prove (fT8T)-(fT9l), it suffices to show RAj — > in part (a), and that 
7 S/=i — > in part (b). To accomplish this, first compute 

RAj — -R5 j^—ipxXxfjTj f 

= sign(aj)S j:R ^_^ x *ti)Tjf, (20) 

where it is permissible here to pass the Riesz transform through the synthesis oper- 
ator because the series for Sj <ih—ti>*\*uTj f converges absolutely in H 1 : the coeffi- 
cient sequence Tjf belongs to I 1 by Proposition[T6l and the synthesizer ip—ip*\*n 
belongs to H 1 by Lemma [191 (noting (fT71) implies © by (H §3.4]). 

Next notice f Rd R(ip — ip * A * /i) dx = 0, since all Hardy space functions and 
their Riesz transforms integrate to zero; cf. (01). Thus in part (b) of the theorem we 
deduce that ± J2j=i RA j ~> °> from <H2> and Theorem CLUb) (and also splitting 

the sum Ylj=i mto two pi eces > where aj > and aj < respectively). 

In part (a) of the theorem we deduce that RAj —> 0, by using (|20l and Theo- 
rem [l7t a), and the following observation. If Pip = 1 a.e., then by computing the 
Fourier coefficients of Pip we find ip(lb~ l ) = for all I G Z d \ {0}, and thus 

(R(iP - iP * X * fi^iib- 1 ) = 0. (21) 

Of course (|2TI ) holds for ^ = too, as observed in the preceding paragraph. Hence 
PR(tp — tp * A * fi) = a.e., by computing the Fourier coefficients of this periodic 
function. 

This finishes the proof. □ 

9.4. Proof of Theorem|8] — synthesis onto H 1 . 

We can assume the dilations expand exponentially, like we did in the proof of 
Theorem [2] 

Let / G H 1 and J G N, and define the sequence cj by (fl~2l where now <ft is 
assumed to be a Schwartz function with <p supported in C$b~ l and f Rd (pdx = 1. 
Then cj G by applying Proposition l2~T1for each j = 1, . . . , J, giving 

1 J 

IMI^fci) < 7 ^ ll^l^llloo||^||i||/||fli = IWIIIoolHiimitfi. (22) 

Observe Scj = ± ^/=i Sj T j/ ^ / in as J ^ oo by Theorem|22] (This 
is where the hypotheses on ip are needed, and that the dilations expand exponen- 
tially.) Combining this with (1221 . we see the open mapping theorem in Appendix lAl 
implies S is open, and that for each e > there exists c G ^(h 1 ) with Sc = f and 

||c||^i(fti) < ||-P|0|||oo||^||i||/||ffi + £• Choosing C = ||P|0| ||oo||^||i proves the 
theorem. 
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9.5. Proof of Corollary [10] — analysis norm for H 1 . 

Take v as in the definition of h 1 in Section 01 Then by boundedness of the 
analysis operator in Proposition ED 

sup H^/IUi < ||pM||oo||^||i||/||hi. 

j 

To prove the other direction of the equivalence, choose ip to be a Schwartz func- 
tion with Pip = 1. Then SjTjf — > / in i? 1 as j — > oo by Theorem [221 a). noting 
that hypothesis (fT7]) is known to hold for the Schwartz function ij) (cf. (6] §3.3]). 
Therefore it follows from Theorem [7] (bounded synthesis) that 

H/llffi <mp\\SjTjf\\ H i < CmvWTjfWf,!, 

j j 

where C is independent of /. 

10. Sobolev space proofs 

Throughout this section we assume the dilations are isotropic and expanding, 
meaning aj = ccjl for some nonzero real numbers aij with |ay| — ► oo. 

We introduce the notation g* r for the convolution of a function g with itself r 
times (for example, g* 2 = g * g), and we write 

A y g = g-g(--y), y G R d , 

for the backwards difference of g by y. Also we define 

I3i(x) = £Oi)l[ ,i)(a?2) • • • l[o,i)( x ^) 

and similarly define functions fo, . . . , (3d, so that the partial derivatives of the unit 
cube indicator function are 

D t t c = A et p t , t = l,...,d. 

10.1. Proof of Theorem [UJ— synthesis £ 1 (u; m «P J a) -> Vy m - P . 

First assume b = I, so that /3 = J<?. Fix c G £ 1 (w m,p , a). 

The initial task is to show that Sc G W m,p . By differentiating the definition ([9]) 
of i/j we convert derivatives to differences: 

D p ij) = {Dxlcy 1 *■■■* (D d l c )* pd * p* m -^ * r? 

= Ag ■ ■ • (/3 X * P1 * • • • * /3* pd * (3* m ~ M *r])£L p . (23) 

Obviously (l23l implies that -0 G W m,p . Further, since P|ry| G L p loc by hypothe- 
sis, we have 

F|/3* w * • • • * /?7 d * /T m ~ |p| * r)\ < P*/ 1 * • • • * (3* pd * /r m ~ |p| * (24) 

and therefore P\D p ip\ G L^ oc by (l23l . This allows us to use as a synthesizer 
when applying Theorem [T] below. 

It is straightforward to show that Sc (which belongs to LP by Theorem [TJ) has 
weak derivatives given by "differentiating through the sum", namely 

D p {Sipc)=S D ^{a\ p \c). (25) 
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Note that the righthand side belongs to L p by Theorem [Q since c G £ (w m,p , a) 
ensures ct' p l c € Thus the function Sc belongs to W m,p , completing the first 

task in the proof. 

The next task is to prove that 

derivatives commute with synthesis, 

in the sense that 

DPS^c = S rip (a^AP C ), (26) 
where we have introduced the function 

(so that for example, ryo = tp). Indeed 

(D p S^c)(x) = E E afc jM (DP^) hk (x) by © 
i>o kez d 

= E a \ Pl E I det °i 1 1/p ( A £ • • • A "» - fc ) b y <E3> 

i>o fcez d 

= E "I"' E ( APc M det aj\ 1/P V P (a>jX ~ k) (27) 
i>o kez d 

by summation by parts, the key step in the proof, 

= \(al"lA"c)(i), 

which proves (l26l) . 

We will deduce an estimate on the synthesis operator of the form 

||5c||i4/™.p < ||-PM||jy(e)||c||4i( l0 m,p )a ), (28) 
which completes the proof of the theorem when b = I. Start by observing 

\\Sc\\ wm .P = W DP Sc\\ P 

\p\<ra 

= E ll^ P (« lpl A P c)ll P by® 

\p\<ra 

< E ii p wiUn^ii alPlA/,c iiw 

\p\<ra 

by Theorem Q] (boundedness of synthesis into L p ). To complete the proof of (l28l . 
notice 

\\P\v p \\\lp(c) < \\P\v\\\lp(c) 

by (HQ). 

To prove the theorem when b ^ I, first rescale the definition © of ip to obtain 
that 

m factors 



M b tp = l c * ■ ■ ■ * tc * M h r\. 
That is, (O holds with M^ip, M^-q and / instead of ift, rj and b, so that M^ip and 
Mbi] satisfy the hypotheses of the theorem for "b = I". 
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By the "6 = I" case of the theorem already proved, then, we have M b ifi G W m ' p 
(so that ip G W m ' p ), and for each sequence c G £ 1 (ty m,p ,a) we have SM b ipjc £ 
W m,p with norm estimate 

ll»S'M i) i/',/ C llvK m .P < ll-P/l-^^lllLP^IIcll^u.^P.a)- 

Further, 

(SM b i>,lc)(x) =^2^2 Cj,k\ detaj\ 1/p (M b Tp)(ajX - k) 
j>o k& d 

= Cj.fcl det Ojl 1 ^! det 6|^(&(ajX — A:)) 

i>o fcgz d 

= | det 6| ^ c^fcl det Oj| 1 / 3 V(a 3 -&a; — 6A;) 
i>o fcez d 

(noting 6 commutes with a,- = ajl) 
= (M b S^ b c)(x). (29) 

Hence 

H^fcSV' 1 & c llw /m ' p = \\SM b ij),ic\\w m 'P 
< \\Pl\M b r]\\\ 

= |det6| _1/p ||P|r/||| LP ( bC )||c|| £ i( w m,p iQ ), 

which finishes the proof. 

10.2. Analysis W m > p w m < p . 

The proof of Theorem [l2]relies on boundedness of analyzers acting on Sobolev 
space, as developed in the next proposition. For simplicity we assume b = I, so 
that the analysis operator at scale j (defined in Section© is just (Tjf)k = (/, 4>j,k)- 

Proposition 23 (Analysis into w m,p ). Assume 1 < p < oo and m G N. Take 
<fi G L q with P\(f)\ G L°°. Fix j > 0, and assume b = I. 

Then Tj : W m ' p — > w m ' p , with norm controlled by the estimate 

J2 {a^WA^fy < \\P\<t>\\U\f\\w^, / G W m > p . 

\p\ < m 

Proof of Proposition \23] Let / G W m ' p and observe 

{Tjf)k = (f, 4> h k) = [ f(y + aj l k)\ det a^'^y) dy. 
Hence for each multiindex p of order < m we have 

al p| A p (T,/) fe = / a f A" /(y + cC 1 *;)! det a^^J^y) dy 

JM. d 

= [ (A p f)(y + aj 1 k) \ det aj\ 1 / q ~cp(ajy) dy 
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where the function A£./ in this last line denotes the p-th backwards difference 
quotient of / with step size aj 1 . Changing variable with y i— ► y — aj k gives 

af^(T 3 f) k = {KJAj,k) = (2}A£./) fc , (30) 

which says that 

differences commute with the analysis operator. 
Thus for each fixed j, taking the £ p -norm with respect to k in (l30l implies 

WalfAPTrfy < ||P|0||| O o||A£./|| p by Proposition [H 

<\\pmu\Dpf\\ p , 

by using the fundamental theorem of calculus. Summing over \p\ < m now proves 
the proposition. □ 

10.3. Scale-averaged approximation in W m ' p . 

Here we prove scale-averaged convergence in W m,p , which we need for the 
proof of Theorem [12] Just like for LP and the Hardy space, the idea is to analyze 
with 4>, then synthesize with tp, and then average over all dilation scales. 

Theorem 24. Assume 1 < p < oo and rj G LP with P\f]\ £ L p oc and J Kd rj dx = 1. 
Lef to € N arcJ define tp by ©. Tafe <ft <E L q with P\<f>\ £ awe? Ld (j)dx = 1. 
Assume b = L. Let f G W m ' p . 
(a) [Constant periodization] If Pi] = 1 a.e. then 

SjTjf -» / in VF m ' p a5 j -» oo. 

[Scale-averaged approximation] dilations aj expand exponentially, then 

1 J 

- SjTj/ -» / in W m ' p as J -> oo. 

The theorem was proved already in our paper [7, Theorem 1], by comparing 
SjTjf with an approximate identity formula. The proof we give below is consider- 
ably easier, and is based on passing derivatives and differences through the analysis 
and synthesis operators. On the other hand, our assumptions on the synthesizer ip 
are noticeably stronger than the Strang-Fix type assumptions in [7], because here 
we assume tp has the special convolution form ©. This is explained in more detail 
in Q Notes on Theorem 1]. 

Theorem [241 a) is due to Di Guglielmo |[28l Theoreme 2'], when rj has compact 
support. Further in this direction, Strang-Fix theory establishes approximation 
rates of the firm 0(\aj\ k ~ m ) in W k,p , for k < m, which improves on the rate o(l) 
in Theorem 1241 a). See Jia [30, Theorem 3.1], and our discussion of the literature 
in (71 §3.5] (where some results special to p = 2 are cited also). 

Theorem l24T b) is needed below when proving Theorem [T2l 
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Proof of Theorem |2?1 Let p be a multiindex of order < m. Then 

D'SjTjf = S^ia^APT.f) by 

•V,.,;/- a;; ./• by© 

= Sj^T^f - D»f) + Sj tVp TjD p f 
= Aj + £?j , say. 

For we have 

P;|| P < WPlVplhpioWTjiA^f - DPf)\\ eP by Theoremffl 

< WPlVpWlLPioWPMlooWKJ ~ D p f\\ p by Proposition [16] 
— ► as j — > oo. 

For Bj, we first note that if Pr\ = 1 a.e. then Pr\ p = 1 a.e. Thus Theorem [171 
implies in part (a) that Bj — > D p f in L p . In part (b), Theorem[T7]implies that if the 
dilations expand exponentially, then j Ylj=i Bj — > in £ p 5 since J Rd rj p dx = 
f Rd rjdx = 1. This proves the theorem. □ 



10.4. Proof of Theorem[T2]— synthesis onto W m ' p . 

In proving S is surjective, we can assume the dilations expand exponentially, 
like we did in the proof of Theorem [2] 

First we prove surjectivity assuming b = I. Let / G W m,p and J G N. Like in 
the proof of Theorem [2] (but with b = I), we take <f> = tc and define cj by (TT2l . 
Then 

ll c j||* 1 (ic> m >P,a!) - (31) 

by using Proposition [23] for each j = 1, . . . , J. And S^jcj —> f in W m,p as 
J — > oo, by Theorem 1241b). 

Hence the hypotheses of the open mapping theorem in Appendix [A] are satisfied 
with A = 1, by (f3TT >. Therefore S^j : ^ (w m ' p , a) —> W m ' p is open and surjective, 
and for each / G l^ m 'P and e > there exists c G f L (io m ' ? ', a) with S^/c = / 
and ||c||^i( w m, P)C j) < ||/||vK m -p + £• This proves the theorem when b = I. 

For the general case where b / /, we rescale like in the proof of Theorem [TT] 
since M^if) and M^r\ satisfy formula ([9]) with "6 = I", the case of the theo- 
rem already proved tells us that for each / G W m,p and e > 0, a sequence 
c G i l {w m,p ,at) exists such that HcH^u^m,}, a ) < ||/||w»™.:p + e and / = Sm^jc 
The calculation (|29T > now implies Mr f = S^^c, as desired. 



10.5. Proof of CorollaryHH — analysis norm for W m,p . 
By boundedness of the analysis operator in Proposition [23] 

sup £ la^HA^/ll^ < ||P|0||U|/||w-m>- 
^ |p|< m 
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On the other hand, choosing 77 as in Theorem l2"4"T a). we see that SjTjf — > / in 
W m,p as j — » 00. It then follows from Theorem QT] (bounded synthesis) that 

< SUP \\SjTjf\\wm, P 

j 

<\\P\ri\\\is<C)Bup E I^I |P| II AP ^/II^- 
7 
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Appendix A. The open mapping theorem 

The open mapping theorem in the following form is used to prove surjectivity 
of the synthesis operator, at various points in the paper. 

Theorem 25. Let X and Y be Banach spaces, and suppose S : X — > Y is bounded 
and linear. Assume 

S(B X {A)) d B Y {1) 
for some A > 0. That is, assume for each y G Y that a sequence {xj} C X exists 
with Sxj — > y as J — > 00 and \\xj\\x < j 4||y||y/ or a ^ J- 

Then S is an open mapping, and S(X) = Y. Indeed, given y G Y and e > 
there exists x G X with Sx = y and \\x\\x < A||y||y + e. 

For a proof, see ll33l Theorem 4.13] with A = 1/5. 

The hypothesis in Theorem [25] can clearly be weakened, to assume only for 
some dense subset of y-values that a sequence {xj} C X exists with Sxj — > y as 
J —> 00 and < A||y||y for all J. 

Appendix B. Discrete Hardy spaces 

In this appendix we study properties of the discrete Hardy space h 1 , which was 
defined in Section |4] We will show h l is independent of the cut-off function used 
in its definition, and that it coincides (when b = I) with the discrete Hardy space 
i? x (Z d ) studied by Q. Y. Sun [37l and C. Eoff 021 in dimension 1, and later by S. 
Boza and M. J. Carro [210 in all dimensions. 

We will need the following result on Riesz transforms of Schwartz functions, 
which is a special case of E Corollary 2.4]. Recall that Z(x) = C d x/\x\ d+1 for 
x ^ 0, and Z(0) = 0. 

Lemma 26. If 6 is a Schwartz function then 

R9{x) = 9(0)Z(x) + O (^pri) f° r al1 xeM. d \ {0}. 

The next theorem replaces the kernel z defining h 1 with a discretization of the 
Riesz kernel, namely the sequence z b = {Z(bk)} having Mi term z\ = Z(bk). 
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Proposition 27. In the definition of the space h , if we replace the sequence z = 
{z k } by z b = {Z(bk)} then we obtain the same space, with an equivalent norm. 

Proof of Proposition \27\ Let v be a cut-off function as in Section 0] Let p, be the 
Schwartz function with /I = v. Note ^(0) = u(0) = 1. Then 

ceo = -*'|§zt| i/ (0 = = ^(0, e g Co, 

where A' is the singular integral operator with kernel 

K(x) = | det b\Z(bx) = | det b\C d 



d+l ' 



|6x| 

Thus by definition of the sequence z in Section|4j we have z k = ((—k) = K/i(k). 
Next, observe that for each x E R rf , 

Kp(x) = | det b\C d p.v. / — tttth(x - y) dy 
Jm d \by\ d+i 

= C d p.v. |— |— -y))dy 

jR d \y\ + 

= R9(bx), 

where 9(y) = ^(b~ l y) is a Schwartz function with 9(0) ^ 0. In particular 
z k = Kp(k) = R9(bk) = 9(0)Z(bk) + O ( * +1 j by Lemma|26] 

= ?(o)4 + o 



Ifcl 



d+l 



for all k ^ 0. Hence 2 — #(0)z fe G i 1 , so that z and z b define identical h 1 spaces 
with equivalent norms. □ 

The proposition and its proof yield a large class of kernel sequences that generate 
h l , for they show that if p is any Schwartz function with ju(0) ^ 0, then we can 
replace the sequence {z k } by {Kfi(k)} in the definition of h 1 . 

Now we can show independence of h 1 from the cut-off function. 

Corollary 28. The space h 1 does not depend on the cut-off function v used to 
define it, and different cut-off functions produce equivalent norms. 

Proof of Corollary |28] This follows from Proposition |27l because z b does not de- 
pend on v. 

Alternatively, consider two different cut-off functions v\ and v-i, giving rise to 
periodic functions (j and £2 as in Section [4] Then £j — (2 is smooth and compactly 
supported in Co and hence has Fourier coefficients in £ l . Therefore the kernel 
sequences associated with v\ and V2 differ by only an i l sequence, and so they 
define the same h 1 space, with comparable h 1 norms. □ 

Corollary 29. Ifb = I then h 1 = H x (JL d ), with equivalent norms. 
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This corollary simply restates Proposition [27] with b = I, because the dis- 
crete Hardy space H 1 (Z d ) is defined (following 0) by the kernel sequence z 1 = 
{Z{k)}\ in other words, 

H 1 (Z d ) = {se£ 1 -.z^sel 1 } 

with a norm ||s||#i/ Z rf\ = 1 1 ^ 1 1 1 — |— 1 1 * s 1 1 1 . Note that in dimension 1, the sequence 
z 1 = {Ijfc^o}/ 71 "^} i s called the Hilbert sequence and was considered by R. E. 
Edwards and G. I. Gaudry lfl4l . who proved boundedness of s i— ► z 1 * s on l v (7L), 
for 1 < p < oo. 

S. Boza and M. J. Carro [3] proved the space H 1 {'L d ) admits a characterization 
by maximal functions in the sense of Fefferman-Stein [ 16], and an atomic decom- 
position in the sense of Coifman-Weiss [11]. The atomic decomposition in one 
dimension was also stated in fTT). It is an interesting problem to investigate these 
characterizations for our space h 1 when b is not the identity matrix. 

Remark on vanishing means in h . If s € h 1 then X^fcez d s & = 0- P roo f : Writ- 
ing a(0 = E feezd s k e- 2 ^ k we see that a(0C(0 = EkM* * z) k e~ 2 ^ k in 
L 2 (Co). This last series is continuous because s * z G I 1 , and <r(£) is continu- 
ous too. But is not continuous at the origin and so <r(0) must equal zero, as 
claimed. 

Conversely if b = I and s E I 1 is finitely supported with X^ez d s fc = 0> men 
s € /i 1 ; cf . (3] Theorem 3.3]. In other words, atoms belong to H. x {JL d ). 

We end this appendix with a question: is h 1 independent of the choice of "trans- 
lation" matrix 6? We suspect not. Of course there is a trivial result: one can always 
replace b by a multiple of b without affecting the resulting space h 1 . 

Appendix C. Banach frames 

This appendix explains how Banach frames arise from the analysis norms earlier 
in the paper. 

Let Y be a Banach space, and let Z be a Banach space whose elements are 
complex sequences indexed by a countable set J. Let {gi}i£i be a subset of Y*, 
the dual space of Y, and let S* : Z — ► Y be a bounded linear operator. We say 
that ({ft}, £**) is a Banach frame for F with respect to Z if the following three 
conditions are satisfied: 

(i) {(/, ft )}eZ,foraU/GF, 

(ii) ||{(/,ft>}||z,faraU/ey, 
(hi) 5*({(/, 5i )}) =/forall/er. 

In other words, "analyzing" with the {gj} maps Y to Z with comparable norms, 
and then "synthesizing" with 5* recovers the identity map on Y. The above defi- 
nition is due to K. Grochenig 11261 : see the treatment in (9j §17.3]. 

The next result reformulates our L p -analysis norm in Corollary [6] as a Banach 
frame result. 
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Corollary 30 (Banach frame for L p ). Assume 1 < p < oo and let <fi and ip satisfy 
the assumptions of Theorem \1 7\with 7=1 (scale-averaged approximation in LP). 
Assume the dilations expand exponentially. 

Then ({| det b\(f)j k}, S*) is a Banach frame for LP with respect to 

1 J 

Z = {c G £°°(£ p ) : S*c = lim — c j,k' , l } j,k exists in LP}. 

3=1 k&Z d 

One can similarly reformulate the Hardy and Sobolev space results (Corollar- 
ies [10] and [15]), for isotropic dilations a 3 - = a>jl. 

Proof of Corollary\30i Write cj = {cj t k}kez d ' so tnat Sj c j = Y.kez d c jP^j,k- 
Define Be = {5jCj}j>o- Then B is a bounded linear operator from £°°(£ p ) 
into £°°(L P ) (by the proof of Theorem [JJ. Therefore Lemma |3~T1 below tells us 
Z is a closed subspace of £°°(£ p ), and hence is a Banach space itself. More- 
over, 5* is bounded from Z into LP . Corollary [6] and Theorem [17] now show that 
({| det b\<pj k}, 5*) is a Banach frame for LP with respect to Z, noting in particular 
by Theorem [171 that 

1 3 

S*(Tf) = lim - V SjTjf = f in L p , for all / G L p . 



*oo J 

3 



□ 



The final lemma states that the preimage of the space of Cesaro-convergent se- 
quences in a Banach space forms a closed subspace. Let Y be a Banach space, and 
write £°°(Y) for the Banach space of all sequences y = {yj}j>o, Vj € Y, such that 

= su Pj>o WyjlW < oo. 

Lemma 31. Let X and Y be Banach spaces and B : X — > £°°(Y) be a bounded 
linear operator. Then the subspace 

1 J 

{x £ X : the limit lim — ^^{Bx)j exists in Y} 



J— >oo «7 



closed in X. 

We omit the proof. 
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